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ABSTRACT 

The time del ay in gravitatio nal lenses can be used to derive the Hubble constant in a relatively 



simple way ( |Refsdal 1964 ). The results of this method are less dependent on astrophysical 
assumptions than in many other methods. For systems with accurately measured positions 
and time delays, the most important uncertainty is related to the mass model used. Simple 
parametric models like isothermal ellipsoidal mass distributions seem to provide consistent 
results with a reasonably small scatter when applied to several lens systems (K oopmans & 
Fassnacht 1999). We discuss a family of models with a separable radial power-law and an ar- 
bitrary angular dependence for the potential ip = r" F(6). Isothermal potentials are a special 
case of these models with (3=1. An additional external shear is used to take into account 
perturbations from other galaxies. Using a simple linear formalism for quadruple lenses, we 
can derive Hq as a function of the observables and the shear. If the latter is fixed, the result 
depends on the assumed power-law exponent according to Hq oc (2 — (3)/ (3. The effect of 
external shear is quantified by introducing a 'critical shear' -f c as a measure for the amount of 
shear that changes the result significantly. The analysis shows, that in the general case Ho and 
7 C do not depend on the position of the lens galaxy. Spherical lens models with images close 
to the Einstein radius with fitted external shear differ by a factor of (3 /2 from shearless mod- 
els, leading to Hq oc 2 — (3 in this case. We discuss these results and compare with numerical 
models for a number of real lens systems. 

Key words: gravitational lensing - distance scale - quasars: individual: Q2237+0305, 
PG 1 1 1 5+080, RX J09 1 1 .4+055 1 , B 1 608+656 



1 INTRODUCTION 

Gravitational lenses offer a unique possibility to de termine cosmo - 
logical distances and hence the Hubble constant ( Refsdal 1964 ). 
The method avoids using any form of distance ladder and is almost 
independent of a deeper understanding of astrophysical processes. 
This has the advantage that possible statistical and systematic un- 
certainties can be controlled much better than with other methods. 
Besides the measurement errors of redshifts, positions and time de- 
lays, the most important source of uncertainties in using gravita- 
tional lenses to determine Ho is the mass model of the lens. Con- 
straints provided by observatio nal data are never suffic ient to fix 
the mass distribution uni quely dSaha & Williams 1997t Williams 
& Saha 2000). Simple parametric models, based on the knowledge 
about typical galaxies in the local universe or desired mathematical 
properties, are normally used to overcome this difficulty. This ap- 
proach, while leading to consistent results, has the disadvantage of 
hiding the underlying uncertainties and making it difficult to quan- 
tify them. 
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In this paper, we consider a family of mass models with sepa- 
rable radial and angular dependence of the potential. External per- 
turbations are included as an external shear. In this way, the dif- 
ferent parts can be studied independently. For the radial depen- 
dence, we choose a simple power-law and are especially interested 
i n the influence the radial slope h a s on the results. Several authors 



dWambsganfi & Paczyiiski 1994 |Witt, Mao & Schechter 1995 



Wucknitz & Refsdal 2001) have found approximate scaling rela- 
tions of Ho oc 2 — f3 for spherical mo dels with external shear. The 
same dependence was found before by Chang & Refsdal (1976) for 



doubles in spherical models without shear (sec also R efsdal & Sur- 
dej 1994). The work of |Williams & Saha (2000| ) showed, that the 
radial slope also has important effects in non-parametric models. 
Due to the very general nature of these models, no exact scaling 
law could be found for them. 

We follow an intermediate approach by using an arbitrary 
function only for the angular part. The models include elliptical 
mass distributions and potentials as well as other shapes. It has been 
observed in several lens systems, that the external shear required to 
fit the data with simple elliptical models is much higher than ex- 
pected. This might indicate a kind of asymmetry of the galaxy that 
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cannot be accounted for by elliptical models. The general angular 
part of the potential we use can describe this 'internal shear' more 
accurately than simple parametric approaches. 

We extend and generalize the work of W itt, Mao & Keeton 
(2000) by including the three independent time delays in quadruple 
systems as constraints for the models. We do not use magnifica- 
tion ratios for several reasons. One is the observational problem of 
reliably measuring correct flux ratios. Fluxes are influenced by mi- 
crolensing and extinction. These effects can be very strong in the 
optical and in some cases they are still significant for radio wave- 
lengths. The other reasons are related to our formalism and will be 
discussed later. 

With this approach, we can derive explicit solutions for the 
Hubble constant as a function of the observables, the power-law 
exponent (3, and the external shear. We use the results to find 
a simple but rigorous scaling law describing the dependence of 
HoAt on f3 in lenses with quadruple images. This scaling law of 
Ho At oc (2— /9)//3 is also valid for more special parametric power- 
law models within the allowed range of /3 and is therefore inherent 
to these models and not an artefact of the general angular depen- 
dence we discuss. 

The qualitative effect of this scaling law is easily understood 
when comparing it with the 2-/3 scaling discussed before and 
with the mass-sheet degeneracy, which leads to HoAt oc 1 — k. 
In all cases, a shallower density profile (larger /3 or k) leads to 
smaller values of HoAt. The flat part of the density distribution 
(k in the mass-sheet degeneracy) amplifies the deflection angles 
but leaves the time delays unaffected. To fit the given geometry, the 
total deflection angles have to be constant, therefore the time delays 
(Ho At) decrease. 

We also discuss the effect of the external shear on the time de- 
lays and the measured Hubble constant. If the shear is changed, the 
internal asymmetries of the mass distribution have to be adjusted 
to compensate for the changes. The total effect of these two con- 
tributions will be described by the new concept of a 'critical shear' 
7 C . The measured Hubble constant is a linear function of the exter- 
nal shear and becomes zero for 7 = <y c . We will present a simple 
interpretation of the critical shear in terms of the image geometry. 

The main goal of this work is to investigate the uncertainties in 
measurements of Ho or cosmological parameters from time delays 
in gravitational lenses. The combination of distances that can be 
determined not only depends on the redshifts and Ho but also on 
the other cosmological parameters (Q and A in standard cosmol- 
ogy) and the smoothness of the mean density distribution ('lumpi- 
ness' parameter a). The results will also help in minimizing pos- 
sible errors, either by selecting lenses with the least uncertainties 
or by using constraints of the model parameters which are lensing- 
independent. 

Finally, we apply the formalism to several known lens sys- 
tems, with and without measured time delays, to compare our ana- 
lytical results with those from numerical model fitting using para- 
metric models. For 2237+0305, we compare with own numerical 
models and find a very good agreement even though the time de- 
lays are not used as constraints in this case. 

Our results can be used directly to determine Ho from time 
delays without explicit modelling but should not be used as a sub- 
stitute for it. They are rather meant to provide an explanation for the 
degeneracies and scaling relations that have been found for several 
families of lens models. Nevertheless, we show that the direct ap- 
plication to real systems is possible. 

In the appendix, we discuss possible Einstein rings in our fam- 
ily of models. Interestingly, infinitely small sources can be mapped 



as elliptical rings for arbitrary values of the external shear. The el- 
lipticity of these rings is directly determined by the shear. 



2 THE LENS MODEL 

For the lens equation and deflecting potential ifi, we use the follow- 
ing notation: 



2 S = Z — Ct(z) 

a(z) = Vtp(z) 

z = r ( cos 9 , sin t 



(1) 
(2) 
(3) 



We use a power-law approach with a general separable angular part 
for the main lensing galaxy. 



^>(z) = r p F{0) 



(4) 



This family of models includes both elliptical mass distributions 
and elliptical potentials with arbitrary radial mass index /3. In the 
following, we assume all image positions as known. For the chosen 
model, this implies also the knowledge of the position of the galaxy 
centre itself, because coordinates are relative to this centre. Later 
we will see, however, that some of the equations are translation 
invariant, leaving the results unchanged when shifting the galaxy. 
The limiting cases of /3 — ► and /3 = 2 correspond to a point 
mass and a density constant in the radial direction, respectively. 
Between these two, the isothermal models are another special case 
with P — 1. The normalized surface mass density of the general 
model is 



a{z) 



f3 2 F{9) +F"(e) 



(5) 
(6) 



Here and in the rest of this paper, primes indicate derivatives with 
respect to 9. A very simple relation holds for the radial derivative of 
the potential which we will need in the time delay equations later: 



z ■ W 



(s) 



(g) 



(7) 



To account for nearby field galaxies or the contribution of 
a cluster, we include an external shear plus an additional con- 
stant mass density or convergence k into our models. This k can 
not be determined from observations as a consequence of the so 
called mass-shee t degeneracy, first discussed by Falco, Go renstein 
& Shapiro (1985) and |Gorenstein, Shapiro & Falco (1988| ). In the 
following, we therefore always use a fixed k. We parametrize the 
shear as follows: 



(7) 



(z) = \z l Vz 



(8) 



The shear matrix T includes both the convergence and the shear 
itself. 



« - 1* ~lv 
-72/ K + 7* 

7 = 7 ( cos 29 1 , sin 29 1 



(9) 
(10) 



With this definition, # 7 points in the direction of the external per- 
turbing mass or equivalently in the opposite direction. Note that this 
part of the potential is a special case of the power-law with /3 = 2. 
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3 TIME DELAYS 

The light travel time t for a certain image at z can be calculated for 
arbitrary lens models by using the equation 



t = -Arff 
c 



D eS = 



D ds 



ip{z)j + const with (11) 
(l + 2d) ■ (12) 



Here Dd and D s are angular size distances of the deflector and the 
source from the observer, while Dds is the distance of the source 
measured from the deflector. We separate this equation as usual into 
one part depending only on cosmology (without Ho), the contribu- 
tion from the Hubble constant, and a lens dependent part. This is 
done by using reduced distances d x (with D x — d x c/ Ho) that only 
depend on the cosmological parameters (and the redshifts) but not 
on Ho itself. We now use a scaled version of the Hubble constant, 

h = ^ , (13) 

which includes the cosmological factors and will later be deter- 
mined by the lensing effect. This h is directly proportional to the 
Hubble constant itself when the cosmology is fixed. For two images 
i and j, we get a time delay of Atij with 

hAtij = h (ti — tj) 

= \ (Kl 2 - KI 2 ) - (V>* - V'j) • 

By using the lens equation (|j|), we can transform this expression 
into a linear functional of the deflecting potential and its deriva- 
tives. Mixed terms like Zi ■ ctj for i ^ j can be eliminated, so that 
the resulting time delay can again be written as the difference of 
the light travel times ti themselves. 



(14) 
(15) 



hti = - i^l + Zi 



V^i -ipi-C 



(16) 



Although it is not immediately noticeable, this equation still is in- 
variant under translations of coordinates. To prove this, one has to 
apply the lens equation. 

For the special potential we discuss here, the relation (^), 
which also holds for ip^ with j3 = 2, can be used to eliminate 
the derivatives of the potential and obtain a simple expression only 
depending on both parts of the potential at the image position: 



-C 



(17) 



We retain using the notion of light travel times for each image in- 
stead of time delays between the images to keep the equations sim- 
ple. The ti are defined except for an overall additive constant that 
gets absorbed into the constant C. 



The last equation was already presented in Witt et al. (200C ) 
for the special cases of (3 = 1 including external shear and for the 
general shearless power-law model. 



Table 1. Parameters for the chosen family of models for a system with n 
images. The Fi are needed in the lens equations. 



parameters 




number 


h 


scaled Hubble constant 


1 


7 


external shear 


2 


P 


power-law exponent 


1 


F 


angular part F(0i) 


n 


F! 


dF/de at 6i 


n 


Z s 


source position 


2 


C 


constant in light travel times 


1 


total 


without fluxes 


2n + 7 



Table 2. Constraints from observations of image positions and time delays 
for a system with n images. 



number 



image positions 
light travel times 



2n 
n 



total 



without fluxes 



3?i 



(uniquely defined) n — 1 time delays without C. Both possibilities 
result in n — 1 more constraints than parameters just for the time 
delays. 

Even for systems with 4 images of one source, it will be im- 
possible to determine all parameters. We therefore decide to fix j3 
in the following calculations so that the results can be used to study 
the dependence of h on the radial mass slope. We will see that, with 
fixed j3, all equations stay linear. 

Another critical parameter is the external shear, which seems 
to be higher than expected in many detailed lens models. Fixing 7 
for the moment, we will be able to investigate the influence of any 
uncertainties in the external shear. In the special case of (3 = 1, it 
will even be possible to determine 7 from the constraints, because 
a number of other parameters do not contribute then. 

We do not include fluxes or magnifications, because they 
would provide only n — 1 constraints (the independent flux ratios) 
and at the same time add n more parameters (the second deriva- 
tives of the angular part of the potential at the image positions F"). 
The effect, that models get less constrained, when more observa- 
tions are included in the analysis, is of course unknown in normal 
parametric models, where the number of parameters is fixed. Our 
models have an infinite number of parameters, of which only a fi- 
nite subset is needed to compare with observations. The number of 
relevant parameters can change when we include more constraints. 

Another reason (besides the difficulties in determining fluxes 
already discussed) is that, in contrast to deflection angles and time 
delays, magnification ratios are non-linear functionals of the lens- 
ing potential, complicating the equations considerably. 



4 COUNTING CONSTRAINTS AND PARAMETERS 

Before solving the equations, we want to discuss how many param- 
eters of the model can at most be determined from observations of 
lens systems with n images, especially n = 4. See Table [I] and ^ 
for a list of parameters and constraining equations. As the time de- 
lays do not change when adding the same constant to all ti, we have 
to include the constant C as parameter and thus have n constraints 
with one parameter for this. We might as well have used only the 



5 LENS EQUATIONS 

To exploit the information given by the image positions, we have to 
insert the power-law potential with shear into the lens equation (JlJ). 
The derivative of the galaxy part of the potential can be obtained 
most simply by rotating its polar form to Cartesian coordinates us- 
ing the transformation matrix 

/a.\ (cose -sme/r\ fd r \ 

\d y J \sind cosO/r ) \dg) ' ( ' 
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Written in a form to take into account the role of 7, Fi and F{ as 
unknowns of the equations, this leads to the following equation. 

-e-i-^G: ;r)(f'M-; 

(19) 

It is easily seen, that this set of 2rc equations can be used to deter- 
mine Fi and F[ , assuming 7 as known: 

PFi = r"' 3 ((1 - «) r 2 - a^ass - 2/i2/ s + 7^ (as? - j/f ) + 2 j y a^J 

(20) 

= r' 13 (yiX B - Xiy B -2j x Xiyi + y y {x\ - ylyj (21) 



6 THE GENERAL SET OF LINEAR EQUATIONS 

We now useJEOj) to express the galaxy potential in the time delay 
equations (|17|), We decide to use the light travel times themselves 
rather than the Ho -independent time delay ratios as constraints. In 
this way we can keep the equations linear and the analysis much 
simpler. As a result, the following set of equations is obtained for 
i — 1 . . .n: 



(1 - K)rl + (x\ - 
2/3 



P 



yi)j x + 2xiyi 7„ 
(Uh + C) = 2]- 



(22) 



P' 



The most interesting fact, besides the linearity, is the simple way 
the mass index (3 appears in the equations. By combining informa- 
tion from the time delay and lens equations the way we did, it was 
possible to remove the terms with /3-dependent exponents. Now (3 
only contributes in the scaling factors of the unknown parameters. 
Having solved the system for one value of (3, we immediately find 
the general solution by scaling h, z a and C with the appropriate 
factors. 



7 THE ISOTHERMAL MODEL 

In the case (3 = 1, the equations ( |22| ) obviously degenerate with 
respect to the source position z s . It is then impossible to constrain 
the latter, but the same information can now be used to determine 
the external shear. In this case, the inclusion of the lens equations to 
determine the parameters Fi of the galaxy potential was not really 
necessary, because f/>' s ' does not contribute to the time delay equa- 
tions (|lj). We can now directly invert the latter or (^) to obtain 
solutions for hi so and 7 iso . 

The transition from exactly isothermal to nearly isothermal 
models deserves some attention. Even for models with (3 differing 
only slightly from unity, 7 is a free parameter, while it is fixed by 
the observational data for (3 = 1. An incorrect external shear for 
almost isothermal models is usually compensated for by diverging 
source positions and Fi and F[, leading to unrealistic models. It 
therefore seems appropriate to use the correct isothermal shear even 
for models that do not exactly obey (3 = 1. However, one has to 
take into account possible measurement uncertainties that introduce 
errors into 7 iso . Especially the time delays, which all contribute to 
the solution, can introduce significant uncertainties. 



8 SOLUTION FOR THE GENERAL MODEL 

In the general case with (3 7^ 1, ( p2[ ) can be solved directly to 
determine h for a given shear 7. Even without writing the explicit 
solution, we see, that 



h oc 



P 



P 



(1 + 5* Ix+gyjy) 



(23) 



for some constants g x and g y . The Hubble constant scales with 
(2— P)/P and is a linear (but not proportional) function of the shear. 
Notably, the isothermal model plays no special role in this equation, 
despite the fact that, strictly speaking, 7 can not be chosen freely 
in this case. As the isothermal shear is usually only weakly con- 
strained due to the limited measurement accuracy, we may, how- 
ever, use different values of 7 even for (3 = 1. This is always true 
when the constraints we used here are not all available. Consider- 
ing this, one might well apply equation ( ^3| ) regardless of the value 
of/3. 



9 THE 'CRITICAL SHEAR' 

As the Hubble constant is linear in 7, there has to be a one- 
dimensional family of values of the shear with vanishing h. The 
shear with the smallest absolute value from this family will now be 
called the 'critical shear' 7 C . If we denote the shearless value with 
ho, we can write the effect of external shear as 



h_ 

ho 



7 -7c 

7c 2 



(24) 



If 7 and 7 C point in the same direction, this scaling factor can be 
written as 1 — j/"/c which is analogous to the scaling of 1 — a/a c 
for an additional mass sheet a with critical density a c . The critical 
shear does not dependent on (3 or the time delays ti and can be 
calculated from the image positions alone. 

If an upper bound of 7 max can be assumed, this translates to a 
range of 



h, 



max \ 
min / 



f ± 



7n 



7c 



/'(> 



(25) 



The critical shear is thus a measure for the amount of external shear 
that can change h significantly. The larger it is, the smaller the 
influence of uncertainties in the shear on the determined Hubble 
constant. Shear in the direction of 7 C contributes maximally, in a 
perpendicular direction the influence vanishes. With the 'direction 
of shear', we mean in this context the orientation 28-, of the vec- 
tor 7. This is not the direction of 6- t towards the perturbing mass, 
cf. eq. (|lo|). External field galaxies located perpendicular to the di- 
rection of maximal effect change h by the same amount but with 
opposite sign. There are, however, four directions where external 
masses do not contribute. See section [lj] for actual numbers of the 
critical shear in real observed lens systems. 



Geometrical interpretation 

For constant light travel times or h = 0, equation (^) describes an 
ellipse whose axes a and b are related to the external shear by 



1-K 



a 2 + b 2 



(26) 



For 7 > 1 — k, this becomes a hyperbola. The position angle of the 
minor axis is the same as that of the perturbing mass responsible 
for the shear (# 7 ). This means, that each ellipse/hyperbola passing 
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Figure 1. Example critical shear -y c (arrow) and values of <y for which 
h/ho = (straight line). Note that the direction of 7 is 20 7 and not the 
direction to the perturbing mass # 7 . 



10 SHIFTING THE LENSING GALAXY 

Surprisingly, most of the quantities we determined do not depend 
on the position of the lens centre. A shift of the centre is equivalent 
to adding a constant displacement to the vectors Zi and z s . If we 
look at the general set of equations (^), we see that such a shift 
adds terms linear in Zi and a constant to the equations. The constant 
term can be absorbed in C and (for 0^1) the linear terms in z s . 
As C and z s are of no interest, the equations do not change when 
applying this shift. 

This means, that h and the critical shear are translation invari- 
ant and can be determined even if the lens position is not known. 
This is only true for the family of models we analyse here. Simple 
parametric models usually only fit the data for a specific position 
of the lens centre. 



-0.5 -1 
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Figure 2. Ellipses passing through all images (stars) of four lens systems. 
The 'roundest' one, which corresponds to the critical shear, is shown as 
solid line. The centres are shown as crosses. For eachellipse, a shear for 
which h/ho = can be calculated by using equation ( 261). a) Q2237+0305, 
b) PG 1115+080, c) RX J091 1.4+0551, d) B1608+656 



through all images cor resp onds to a value of 7 with h/ho 
According to equation (E4), for these values 

7 ■ 7c = 7c 



(27) 



holds. The values of 7 for all these conic sections span the complete 
subspace for which h/ho = (see Fig. The critical shear is 
the smallest of these values and can therefore be calculated with 
equation (26J) for the 'roundest' ellipse passing through all images. 
In Figure Iwe show several such ellipses for some lens systems. 

We conclude, that for certain values of 7 (see Figure |l|), e.g. 
for 7 = 7 C , all time delays vanish. Each of the fitting ellipses can 
act as an Einstein ring for the correct value of 7. The light travel 
time is the same for all parts of such rings as Fermat's theorem 
requires. In our consideration, the potential is only constrained at 
the image positions and arbitrary in other directions. The Einstein 
ring may therefore break up and form a number of discrete images 
with still vanishing time delays (see appendix). 



11 SPHERICAL MODELS FOR NEARLY EINSTEIN 
RING SYSTEMS 

For spherical models, the equations are overdetermined. There are 
nevertheless systems, which can be fit accurately with this kind of 
model. It can be shown, that h sp h/ho becomes zero for — > to 
cancel the vanishing denominator in ( p3| ) and ensure a finite h Bp h. 
It can also be shown, that 7 is parallel to 7 C for arbitrary 0. Taken 
together, this means that for point mass models, we always have 
7 = 7 C . The geometrical interpretation can be used to determine 
the direction of external shear for spherical models without calcula- 
tions. It is given by the orientation of the minor axis of the roundest 
ellipse passing through all images. For the systems discussed be- 
low (and for point mass lenses), even the absolute value of 7 can 
be determined from this ellipse. 

For systems, where the n images are all located very close to 
the Einstein ring at r , we can recover another scaling relation. In 
this case, the power-law can be interpreted as a local approximation 
to any radial mass profile, like softened power-law spheres or other 
models. We assume, that one fitting reference model is known. It 
is then possible to find a family of other models which are also 
consistent with the observations. 

TJiis_^as_dOTie_jiumerka^ & Paczyhski 

( 1 994) a nd led to a scaling of /t SD h oc 2 — /3. Wucknitz & Refs- 
dal (2001) presented a simple interpretat ion of this fact in term s of 
the well known mass-sheet degeneracy (Gorenstein et al. 198t). If 
we multiply the lens equation z s = z — a with 1 — k, we get an- 
other lens equation with the source position and lensing potential 
(or deflection angle) scaled with the same factor, plus an additional 
constant convergence k. This means, that lens models given by the 
deflection angles a and a = (1 — k) a + nz are equivalent, but 
source positions, time delays, etc. are scaled by 1 — k in the latter 
model. 

As reference model, we choose an exponent 0o = 1. This 
model is now transformed as just described and then approximated 
locally near ro by a modified power law with exponent 0, Figure ^ 
illustrates this idea. For the new exponent, we find 



= 1 



(28) 



As time delays, source position and external shear all scale with 
1 — ft, we find the general scaling laws 



z s oc 2 — 0, h Bph oc 2 — 0, 7 oc 2 — 



(29) 



for an arbitrary mass index 0. The time delay ratios do not change. 
With 7 = 7 C for point mass systems, the shear for arbitrary 
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a = (1 — re) + re r 
& = r 1 / 2 

a = l 



1 1 1 1 

1 2 

r/r 

Figure 3. The deflection angles for the three equivalent lens models. The 
isothermal (f3 = 1) reference model is shown solid, the transformed model 
dashed and the final equivalent power-law model dotted. 



(3 becomes 7 = 7 C (2 — /3)/2, leading to a ratio of the Hubble 
constants for the spherical and shearless case of 

The two models differ significantly for a realistic range of 13. 



12 INFLUENCE OF THE RADIAL MASS INDEX /3 

One of the interesting properties of equation (^) is the very simple 
dependence of the results on (3. The determined value of the Hub- 
ble constant h simply scales with the factor (2 — (3)/f3. The most 
alarming fact is, that this factor does not depend on the geometry 
of the lens, the time delay ratios or the amount of external shear. 
When using the models described here to determine Ho from lens 
time delays, the error due to the assumption of an incorrect f3 will 
be exactly the same for all lens systems as long as the real f3 is more 
or less equal for all le nsing galaxies. 

Witt et al. (2000) numerically found a scaling of (2 — 13)/ (3 in 
the case of a power-law model without external shear for orthogo- 
nal image pairs (#; — Qj ~ 90°), while these computations lead to 
2 — P for opposed images. In contrast to their work, we have used 
all time delays as constraints so that they cannot scale differently. 
The common scale factor of all time delays shows as a scaling of h 
in our calculations. 

The reader might feel as uneasy about the seemingly diverging 
time delays or h in the limit (3 — > as the author did. This limit is 
equivalent to point mass models and one should not observe diverg- 
ing time delays for this kind of lens. The point causing trouble here 
is the fixed external shear in our considerations. To fit the data with 
a point mass model, the shear has to be equal to the critical shear 
7 C . Taking this into account, the result gets multiplied by a vanish- 
ing h/ho oc j3 which cancels the 1/(3 factor. If we now change the 
shear by a small amount, the potential and thus h will immediately 
diverge in the limit of small {3. This has no direct physical impli- 
cations, because the mass models will become extremely unphys- 
ical to compensate for the shear effect. Small relative differences 
in the Fi, that will be introduced by an incorrect shear, will lead to 
enormous asymmetries in the mass distribution. This is related to 
the fact, that realistic compact mass distributions can provide only 
almost spherically symmetric potentials. Any multipole moments 
would radially decrease more rapidly than the monopole term and 
must be very strong to have any ef fect. We will discuss this in more 
detail for a special case in section 13.1 . 



For spherical systems with images near the Einstein ring, we 
confirmed the approximate scaling of h sp h oc 2 — (3. This seems to 
be incompatible with the general scaling law of ho oc (2 — (3)/ (3 at 
first sight. With the factor of h Bp h/ho = (3/2 from equation (f30|), 
by which the shear changes the result in spherical models, these 
two results are, however, in perfect agreement. 



13 APPLICATION TO SPECIAL CASES 

To illustrate the results, we have presented here, and to test their 
relevance for real lenses, we want to apply the formalism to sys- 
tems with a special Einstein cross like symmetry and to some of 
the known real systems that are either useful to actually determine 
H (1115+080, 0911+0551, 1608+656) or are interesting because 
they are very well studied systems like 2237+0305. The calcula- 
tions will show, that both scaling relations (/i sp h oc 2 — (3 and 
ho oc (2 — (3)/ (3) are relevant for the determination of Ho from 
time delays. The detailed numerical models for 2237+0305 will 
furthermore show, that the scaling also applies if the time delays 
themselves are not used to constrain the models. All time delays 
scale almost exactly as predicted by our analytical work. We will 
also see that parametric models may fit only for a limited range of 
(3. The scaling relations are then valid only within this range. 

Besides the effects of the exponent /3, the possible strong ef- 
fects of any external shear will also be confirmed by the numerical 
models. 



13.1 Symmetric Einstein cross like systems 

A rather special example of systems shall be discussed explicitly in 
this section. We consider a lens with time delays Ati2 = At$4 — 
and the following image positions: 

xi,2 =±(l-e)ro yi,2 = (31) 
£3,4 = j/3,4 = ±(1 + e) r (32) 

From these data, we immediately conclude z s = 0. The y- 
component of the shear does not contribute at all and cannot be 
constrained. The equation determining the Hubble constant now 
reads 



hAt 13 = (l~K)r 2 ?— P-leU- 1 - 

P V 7c 



with the critical shear 



7c = (1 - k) 



2s 



(33) 



(34) 



1 +e 2 

The special symmetry makes it possible to choose any value of ex- 
ternal shear even for isothermal models. Furthermore, it is possible 
to exactly reproduce the data with spherically symmetric models 
plus external shear. In this case, the shear is uniquely defined: 

7 = (1 - k) (2 - (3) e + 0(e 2 ) (35) 

The time delay equation now becomes somewhat more complicated 
than in the non-spherical case. To first order in e, it reads 



Ati 3 = (1 - it) rl (2 - 0) e + 0{e 2 ) 



(36) 



Comparing this with ([33]), we recover the factor of (3/2 between the 
spherical and shearless case. Both models fit the data exactly and, 
in this special case, the models are even compatible with highly 
popular elliptical mass distributions. That means that without any 
independent information about the external shear (or equivalently 
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the ellipticity of the galaxy itself), we have a factor of two uncer- 
tainty even when only considering these two simplest models for 
(3 = 1. The real situation may be even much worse, when we con- 
sider models with internal and external shear. In this case, any small 
unknown contribution of external shear of the order of -y c (which 
for very symmetric systems becomes arbitrarily small) will change 
the result significantly , 

Witt et al. (tyso) discussed exactly the same type of systems 
with spherical models plus shear. With 7 fixed, they also derived a 
scaling law of (2 — (3) j (3 (see their equation When the shear is 
constrained by the lens equations, the scaling changes to the 2-/3 
form. 

We finally want to discuss the consequences of diverging time 
delays in the fixed shear case for (3 — *• due to equations (E3() and 
(^). For simplicity, we assume 7 = 0, re = and e < 1, but 
the argument is generally also true for other values. We write the 
potential as a multipole expansion)^] 



/3F{8) = rf' 13 ( 1 - ^ a fc cos 2k9 



(37) 



To be compatible with equation (po|), the coefficients have to meet 
the condition 



E 



a k = (2-/3)e 



(38) 



We notice, that densities (g) can become negative near the axes. To 
minimize the angular density contrast, we have to keep only the 
monopole and quadrupole terms and set all higher coefficients to 0. 
The potential is then equivalent to a density of 



1 cos 26 



(39) 



which is everywhere positive only for sufficiently high values of 
(3, When using realistic mass models, we can therefore expect a 
lower bound for /3 to achieve acceptable fits. This applies not only 
to this special symmetric lens system but is true in general. Numer- 
ical models presented in the next section will confirm this result for 
7 = (see Figure m. 



13.2 The Einstein cross Q2237+0305 

This lens is not usually taken into consideration when thinking 
about determination of Ho, because the time delays are expected 
to be very small and can therefore not be determined easily. Here 
we show that even if all three time delays were known exactly, con- 
straints for the Hubble constant would still be very weak. 

The degenerac y caused by the unknown mass index j3 was 
already discussed by WambsganB & Paczynski (1994 ) for spherical 
models plus external shear. The authors found the 2-/3 scaling 
using numerical models. We now want to investigate how strong the 
assumption of a spherical main galaxy really influences the results. 

No time delays are available for 2237+0305 and they may 
never be determined. We can nevertheless calculate the critical 
shear defined before and compare it with the a typical value one 



1 The exponent of the first term in equation (8) in Witt et al. (1995) is 



incorrect, it should be the same as that in the second term. (S. Mao, private 
communication) 

2 We use an expansion for the principal axes of the lens system. The 
sin2fc0 terms might be included as well but they would not change the 
density on these axes, which is what we are interested in. 



gets for spherical m odels. Positions inc luding error bars used for 
this were taken from Crane_££ j aL_(199l|) to make results compara- 
ble with WambsganB & Paczynski (1994). 

The critical shear as calculated from these positions^] is 7 C = 
(0.0915 ± 0.044, 0.0958 ± 0.045) or j7 c j = 0.13 where the er- 
rors are la bounds from Monte Carlo simulations. Numerical mod- 
elling results in a shear of 0.0696 almost exactly parallel to 7 C for 
isothermal spherically symmetric potentials. We therefore expect 
the time delays (or h if we take At as known) of the spherical 
model with shear to be a factor 0.4745 smaller than in the shear- 
less case. For the moderately small shear of w 0.07, this is a huge 
effect. This factor is in good agreement with the expected value of 
1/2 for idealized systems. 

To compare results in the general case, we performed numeri- 
cal model fitting with an elliptical potential approach plus external 
shear. 



F(6) oc 



1+e 



+ 



1-e 



(40) 



Elliptical potentials are known to be unphysical for large elliptic- 
ities. Although e is small in our case, we may expect unrealistic 
solutions for small values of (3, because the limit of acceptable el- 
lipticities vanishes for [3 — > (cf. last section). In fact the fitted e 
also increases with decreasing (3. 

We decided not to use invented time delays (calculated for a 
reference model) to fit the models. In this way, we can check the 
validity of our results even for cases where multiple time delays 
are not used as constraints. Plots of the residuals, ellipticity, shear 
and time delay between component A and B are shown in Figure^. 

The non-vanishing residuals at (3 — -> 2 might be worrying at 
first, because a spherical model without external shear can fit any 
image configuration for (3 = 2. The result would be a sheet of con- 
stant density equal to the critical density. There would be no iso- 
lated images but an area of constant (and very high) surface bright- 
ness. One might thus naively think, that residuals should be very 
small near (3 = 2. This is not the case. It is true, that the deviations 
of the projected images will become arbitrary small in the source 
plane. On the other hand, however, the magnifications diverge in 
the limit, causing the deviations in the lens plane to stay finite. 

Three families of models have to be discussed. First, we fixed 
the ellipticities at 0^o^om£are_jvifh^he^ej^ 
& Paczynski (1994). The residuals are almost constant for spherical 
models. The shear and the time delays scale very accurately with 
2 — (3 as in the idealized considerations. 

More interesting in our context is the behaviour in the shear- 
less case where we expect to find a (2 — (3) / f3 scaling of the time 
delays. The \ 2 S ets unacceptably large for (3 <C 1. In the other 
cases (J3 > 0.6), the agreement with the theoretical predictions 
also shown in Fig. ^ is very good. For isothermal models, the ratio 
of time delays calculated for the two models (spherical/shearless) 
is 10-20 per cent larger than predicted by the critical shear. This is 
still a very good agreement considering that the time delays were 
not used as constraints for the numerical models. 

Real lensing galaxies usually are elliptical and also embedded 
in an external shear field. We therefore also let both e and 7 vary 
freely in order to minimize % 2 ■ Counting the formal number of con- 
straints and parameters, we expect a minimum of % 2 = 0. The fact, 
that x does not vanish is a confirmation of th e degeneracy involv - 
ing shear and ellipticity already discussed by Witt & Mao (1997). 



All coordinates in this paper: x to east and y to north. 
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1.5 



ellipt.pot. shearless EP 
ellipt.pot. + ext.shear EP+XS 
spher. + ext. shear XS 




posit ions, does not affec t the result. Using the HST observations 
from Impey et al. (1998) with their claimed accuracy of 0'.'002 as 



Figure 4. Results of numerical model fitting. The top panel shows the re- 
duced x 2 f° r the best models. Here f3 was counted as fixed, resulting in 
v = 3 / 3 / 1 degrees of freedom for the spherical / shearless elliptical po- 
tential/general elliptical potential model, denoted as XS/EP/ EP+XS. In 
the middle, both components of the shear and ellipticity are plotted for the 
EP+XS model. The last panel presents the time delay A— B in hours for 
Ho = 75 kms" 1 Mpc" 1 and Einstein-de Sitter cosmology. We also in- 
cluded a curve scaling with (2 — /3)//3 to compare with the EP case. 



The effective number of parameters is therefore smaller then the 
formal one. Because we notice the ellipticity changing only slightly 
with f3, a scaling like that in the spherical case is expected. This 
can indeed be seen in the plot, where the time delay scales with 
2-/3. Contrary to the spherical case, 7 does not scale proportional 
to 2 — P but is additionally shifted by a constant offset. This is due 
to the fact that part of the shear has been transformed to ellipticity. 

We conclude, that both scaling laws for h or the time delays 
can be relevant, depending on the family of models used. In the case 
of 2237+0305, however, the influence of external shear is stronger 
than the effect of (3 for any realistic values of the latter. 



13.3 PG 1115+080 

Models and their degener acies for this quadrup l e system have been 
studied extensively (e.g. |Courbin et al. 1997|; fmpey et al. 1998[ 



Keeton & Kochanek 1997|; Saha & Williams 1997j; Schechter et al. 
1997; [Zhao & Pronk 2001 ), leading to a variety of more or less 
realistic mass distributions and a range of values for the Hubble 
constant between about 40 and 80 kms -1 Mpc -1 . 

All these authors agree on the importance of taking the ef- 
fect of a nearby galaxy group into account. Our formalism in- 
cludes this group as an external shear of the order 7 ~ 0.1. This 
shear is however not well constrained. Keeton, Kochanek & Seljak 
(1997) showed, that the residuals do not change much for ranges 
of 7 ~ 0.06-0.2. With our general family of models, the effect 
of unknown external shear can be quantified by using the critical 
shear and equation (^). To calculate 7 C , we only need image po- 
sitions relative to an arbitrary reference centre. The uncertainty in 
the galaxy position, which is usually much higher than in the image 



a basis for Monte Carlo simulations, we obtain a critical shear of 
7 C = (0.142 ± 0.002, 0.167 ±0.003 ). Although the ground based 
positions from Courbin et al. (1997) are not compatible with the 



HST results within the formal error bars, the critical shear from 
these data is about the same, 7 C = (0.143±0.003, 0.160±0.005). 
This means, that any shear of the order 0.1 can change the results 
for Ho significantly. As the uncertainties in 7 are of this order of 
magnitude, large effects on h result even for fixed (3. 

The formalism we present in this article was not developed 
to directly determine Ho from observations but to study the model 
degeneracies and scaling laws. One might nevertheless try to use 
the prescription from section to obtain an estimate for the Hubble 
constant and the external shear for an isothermal model. The errors 
of the observational data of course have to be taken into account. 

Time de lays derived from the same light curves ha ve been 
published by |Schechter et al. (1997| ) and |Barkana (1997| ). The A 
component was not resolved in these observations, and time delays 
were only determined relative to the sum of Ai and A2. This is jus- 
tified by the small time delay between the two, which is expected 
to be the order of hours. For our Monte Carlo simulations, we as- 
sumed a time delay between the A images of (0 ± 1) days. 

Results for Ho in the isothermal case differ depending on 
which set of time delays and positions i s used. With the red- 
shifts z s = 1.722 and z d = 0.310 |Tonry 1998\ and an 



Einstein-de Sitter universe, we obtain values of Ho between 47 and 
58 kms" 1 Mpc - with errors between 12 and 30 per cent {la). 
The external shear is only weakly constrained, but seems to be of 
the order 0.1. We have to stress, that this result includes all pos- 
sible isothermal models with arbitrary angular dependence and is 
thus much more general than elliptical mass distributions. As long 
as different determinations of the positions and time delays are not 
consistent with each other within their error bars, any results for 
Ho have to be interpreted with care of course. 



13.4 RX J0911.4+0551 



This quad, initially discovered as a triple (Bade et al. 1997), has 
unique geometrical properties and is a strong candidate for time 
delay determination, although no result has been published yet. 
Rapid variability has been detected in the X-ray regime, provid- 
ing the possibility of a d etermination of all t hree time delays with 
unprecedented accuracy ("hartas et al. 2001 ). Our first models (el- 



liptica^jjotottia^hjs^iea^^ et al. 

(1998). The external shear in the best-fitting model is 7 = 0.32 
and points almost exactly in the direction of a nearby cluster of 
galaxies. 

The redshift and velocity dispersio n of this cluster was mea- 
sured by [Kneib, Cohen & Hjorth (200fj|) to z = 0.769 and a v = 



(836 t"20o) kms -1 . From this, they derive an absolute shear of 
7 = 0.11 t-o oi usm S a SIS model for the cluster. 

To compare with lens models, we have to use the reduced 
shear 7 = 7/(1 — «), because the convergence n caused by the 
cluster was not taken into account explicitly. The mass-sheet de- 
generacy simply scales all parameters in (|2^) with 1/(1 — k). For 
SIS models, 7 = k holds and we obtain 7 = 0.12 1" 'q4 . This mea- 
sure ment is not in good agreement with the model from B urud et al. 
(1998). The reduced shear differs by about 0.2, the direction of the 
two being almost identical. A possible explanation for this discrep- 
ancy is the presence of a second galaxy close to the main lens that 
might change the potential considerably. It is also possible that the 
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internal asymmetry of the main galaxy itself can not be described 
as an elliptical mass distribution. 

To estimate the uncertainty in time delays or the Hubble con- 
stant derived from them, we use again the critical shear which is 
7 C = (-0.553 ± 0.013, 0.101 ± 0.005) or j 7 j = 0.56. Even with 
this very large critical shear, the uncertainty in the real 7 has sig- 
nificant effects, because it is large as well. 



13.5 B1608+656 

This system is the first and up to now only quad for which all 
three jndejsend^nttmiedekr^^ et al. 

1999). This offers the unique possibility to apply our method to 
a system providing the complete set of constraints. HST images 
show a main lensing ga laxy but also a wea ker second galaxy be- 
tween the four images (lackson et al. 1998). We nevertheless ap- 



ply the method to B 1608+656, if only to see if the effect of a sec- 
ondary lens can be detected in this way, e.g. by pretending there is 
a very large external shear. Data for the positions of the images and 
the main lensing galaxy were taken from K oopmans & Fassnacht 
(1999). The formal accuracy of the image positions is extremely 
high, of the order 2-12 |i.arcsec. For the Monte Carlo simulations, 
we used 1 mas scatter in each coordinate to account for possible 
shifts by local density fluctuation s, caused for example by globular 
clusters ( jvlao & Schneider 1998 1. 

We can use the general equation (|2^) to determine the Hubble 
constant for shearless models. With the redshifts of Zd = 0.6304 
and z s — 1.394 and standard Einstein-de Sitter cosmology, we 
obtain a value of H = (37 ± 5) kms" 1 Mpc" 1 for (3 = 1. 
The critical shear is f c = ( 0.072 ± 0.001 , 0.069 ± 0.001 ). For 
isothermal models, equation (^) predicts a shear of about 7 = 
( -0.32 ± 0.02 , -0.11 ± 0.01 ), depending on which HST image 
is used to determine the galaxy position. The result for the Hubble 
constant is H = (130 ± 15) kms -1 Mpc -1 . 

The enormous differences in both Ho predictions is a conse- 
quence of the large external shear 7 2> 7c . No external shear at all 
is needed to fit the data when inclu ding the influence of the second 
galaxy in the field. The models in Koopmans & Fassnacht (1999) 



even predict velocity dispersions for both galaxies which are of the 
same order of magnitude. This is surprising, since the secondary 
galaxy is much weaker in all bands in the optical images. 



14 DISCUSSION 

We used a very general semi-parametric lens model approach to 
study the changes of time delays and the determined Hubble con- 
stant with the assumed radial density slope, quantified by the radial 
mass index /3. By using only linear constraints, it was possible to 
keep all the fundamental equations linear. This made the study of 
the p dependence easy. For fixed external shear in quadruple lenses, 
this resulted in the simple scaling law Ho oc (2 — 0)/f3, indepen- 
dent of the lens geometry, the time delay ratios or the external shear. 
This means, that a systematic error in the assumed will have ex- 
actly the same effect on all lenses and will not show as scatter in the 
results. The g ood agreement between Hp me asurements from dif- 
ferent lenses ( Koopmans & Fassnacht 199S| ) should therefore not 
be taken as evidence for an accurate determination of Hq. It merely 
shows, that all lensing galaxies seem to have more or less the same 

a. 

In nearly isothermal models, a systematic error of only 
10 per cent in will result in an error of about 20 per cent in the 



deduced Hubble constant. To compare the results from lensing not 
only with each other but also with results derived from other meth- 
ods, this possible source of error has to be taken into account. 

Furthermore, it is important not only to be aware of this effect, 
but to try and obtain better constraints on the radial mass profile. 
Possible ways to do this include modelling lenses with multiply 
imaged extended sources to obtain constraints for the lensing po- 
tential at a wide range of distances to the lens centre. Lenses with 
multiply imaged point sources can constrain the models only at a 
small number of radii. This should be done not only for systems 
with measured time delays, but for as many applicable lenses as 
possible to acquire reliable measurements of /3 for a representative 
sample of lenses. 

Another possibility is the detailed study of the dynamics of 
lensing galaxies or of other (usually low redshift) galaxies to learn 
more about the range of realistic galaxy mass profiles and use the 
results of this analysis to produce more realistic models for the 
lenses. 

We also quantified the effect of external shear by introducing 
the concept of a 'critical shear' 7c . The effect of 7 on Ho is lin- 
ear and strongest in the direction of ~f . For a fixed direction, its 
amount is proportional to 7/70- The shear has to be exactly crit- 
ical to fit the observations in point mass lenses. The value of 7 C 
can be found in a geometrical way. It is given by the ellipticity of 
the roundest ellipse passing through all images. For 7 = 7 C and 
non-vanishing /3, the time delays HoAt become zero. This is also 
true for a whole family of models which are represented by the less 
symmetric ellipses fitting the images. In the appendix, we discuss 
the relation to possible Einstein rings which are exactly given by 
these ellipses. 

The effect of shear is also the clue in understanding the com- 
patibility of the general scaling law with the simpler one of ft sp h oc 
2-/3 for spherical models. In the latter models, the shear is con- 
strained by the observational data and changes ft by a factor of 
[3/2 when compared to shearless models. In cases where spherical 
models are able to fit the data, the allowed range of 7 results in an 
uncertainty of h always covering both models. This may in certain 
cases only apply for a limited range of /3. 

Interestingly, the value of the critical shear and the Hubble 
constant in the general model (with fixed shear for Ho) do not de- 
pend on the position of the lensing galaxy. This may be of use for 
systems where this position cannot be determined accurately. 

When using more general models or less constraints than in 
our calculations, the scaling laws still apply. They are only valid 
for a subset of the possible models then and one would have to 
expect even larger uncertainties when using the whole set. This also 
applies for lenses with less than four images. More special models 
on the other hand, like parametrized elliptical power-law models, 
may be able to constrain the range of possible results much better. 
Nevertheless, the scaling laws still apply for the range of models 
that are compatible with the constraints. Even in these cases, our 
results may be used to determine Ho without explicit modelling. 
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APPENDIX A: EINSTEIN RINGS AND HIGH IMAGE 
MULTIPLICITIES 

An interesting property of the general power-law models we used 
in the main part of this paper is the possibility to produce Einstein 
rings from point sources for arbitrary values of the external shear. 
For an Einstein ring parametrized by r = r(9), all points on this 
ring must have the same light travel time to meet Fermat's theorem. 
For this appendix, we set k and 7„ to for simplicity. For an arbi- 
trary direction of the shear, we just have to replace 8 by 8 — # 7 . We 
also assume, that 7 < 1. 

With t = 0, the general equations (^) describe an ellipse, 
which is centred on the lens in the special case z& = or /3 = 1. 



The minor and major axes are ro / ■s/l ± 7, compatible with equa- 
tion (|2q). The potential is in this case an elliptical one 



F(8) = i^" ,3 (l+7 cos26») /5/2 

P 



(A2) 



With respect to the tangential caustic, the effects of ellipticity 
and shear cancel in these models and the caustic degenerates to 
a point. This is qualitatively different from elliptical mass distribu- 
tions where the caustic is deformed and overlaps itself, producing 
areas of higher multiplicities (see below). With arbitrary z a and /3, 
the centre of the ellipse is shifted to zq witbjj] 

1 + 7* Iv \ r _ 1-/3 
7v 1-7*/ 2-/3 



(A3) 



For P ^ 1, this shift can take any value if z s is varied. To obtain 
a globally unique function F(6), the centre of the lens has to be 
located inside of the ellipse. 

Even for lens systems with four images, it is always possible 
to find an ellipse passing through all of them, which can act as an 
Einstein ring for the corresponding value of 7 given by the ellip- 
ticity. This does not mean, that we always see an Einstein ring for 
this special value of external shear, as F(8) is not constrained for 
angles between the images. 

Small deviations from the Einstein ring case can lead to an 
arbitrary number of images near the former elliptical ring. Special 
cases of these systems (singular isothermal ellipsoidal mass dis- 
tributions with shear) with up to eig ht im ages were mentione d by 
Lopez, Wucknitz & Wisotzki ( 1998|) and |Witt & Mao (200C| ) and 
discu ssed in detail by |Keeton, MacT& Witt""(2000|). Evans & Witt 
(2001) present results for shearless models with arbitrary F(9). 
From (p()) and ( |2l| ) we obtain the following condition for an im- 
age of a source at z s — 0: 



F[_ 
F 



pj sin 28 
1 + 7 cos 28 



(A4) 



A global solution for this differential equation is given by (A2) 
which leads to the elliptical ring we discussed before. For a number 
of discrete images, a more general solution is possible: 



F{8) = f(8) (1+7 COS260" 72 



(A5) 



At the positions of the images 9i, ( A4 ) has to be met, leading to the 
simple condition 



(A6) 



As f(9) is an arbitrary function, we can easily construct systems 
with any number of images. The radial coordinates of the images 
can then be determined to be 

r{0i) = (/3/(0 l )) 1/<2 ~' 3) (l + 7Cos20 l )~ 1/2 ■ (A7) 
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author. 



r = r (1 + 7Cos26>)" 1/2 



(Al) 



4 This equation is valid for arbitrary directions of 7. 



